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Abstract 

Recently, a Monte Carlo method has been presented which allows for the form- 
free retrieval of size distributions from isotropic scattering patterns, complete with 
uncertainty estimates linked to the data quality. Here, we present an adaptation to 
this method allowing for the fitting of anisotropic 2D scattering patterns. The model 
consists of a finite number of non-interacting ellipsoids of revolution (but would work 
equally well for cylinders), polydisperse in both dimensions, and takes into account 
disorientation in the plane parallel to the detector plane. The method application 
results in three form-free distributions, two for the ellipsoid dimensions, and one for 
the orientation distribution. It is furthermore shown that a morphological restriction 
is needed to obtain a unique solution. 
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1. Introduction 

Anisotropy in materials nanostructure is often induced upon fabrication or process- 
ing into industrially relevant shapes (e.g. wires, fibres, bars and sheets). For high- 
performance polymer fibres, for example, the oriented structure is effected by fibre 
spinning and further refined in tensioned heat-treatment procedures (Northolt k, 
Sikkema, 1991). Such an oriented structure has been shown to have a high corre- 
lation to the physical properties of the material (Rao et al, 2001). It is therefore of 
significant interest to extract as many structural parameters as feasible to achieve a 
good insight in the structure-property relationships for such anisotropic morphologies. 
Small-angle X-ray scattering can obtain statistically relevant bulk structural param- 
eters for a plethora of such materials. 

In small-angle scattering these anisotropic morphologies can give rise to anisotropic 
scattering patterns, which inherently contain an increased dimensionality of informa- 
tion. When such scattering patterns are sectioned, with each section treated as if 
it originated from an isotropic structure, invaluable cross-correlations are invariably 
lost. Fitting of a complete two-dimensional anisotropic scattering pattern is therefore 
of high interest. Prior work by others (Fischer et al, 2010; De Geuser et al, 2012) 
as well as this author (Pauw et al, 2010) on describing such anisotropic scattering 
patterns based on classical 1 descriptions of data have led to useful but complicated 
methodologies which tend to be rather inflexible and require tuning to the sample at 
hand. Development of a more flexible, generally applicable method is therefore rather 
important. 

While for isotropic scattering patterns several such generally applicable methods 

"Classical" here implies the combination of a form factor with a structure factor and predefined size 
distribution, ofttimes fitted using a least-squares minimisation procedure commonly with up to 10 
variables, as described in detail by (Pedersen, 1997). Conversely, the "generally applicable methods" 
in this paper describe methods that retrieve (form-free) either the form factor, structure factor or the 
size distribution when the remaining two are assumed. 
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have appeared over the years (Glatter, 1977; Magnani et al., 1988; Krauthauser 
et al., 1996), for anisotropic scattering patterns such form- free retrieval methods have 
been much more slow to develop. Of these, the inversion method by Stribeck (2001) 
in particular has been a significant milestone, capable of retrieving two-dimensional 
correlation function for anisotropic scattering patterns. Despite providing information 
in real space, however, correlation functions remain difficult to interpret. 

Shown here, therefore, is a modification to our recent MC method (Pauw et al, 2013) 
to allow for full fitting of anisotropic scattering patterns. This modification leads to 
the extraction of three form- free distributions, with two describing the anisotropic 
scatterers' dimensions and the third the orientation distribution. The first application 
results are discussed after a brief description of the necessary changes implemented 
to its isotropic forerunner. Parameter retrieval from a simulated scattering pattern 
shows that the original parameters are reasonably well retrieved. We think these are 
encouraging results and by presenting them we hope to foster discussions towards 
further development of the proposed MC method. 

2. Required adaptations to the MC method 

The MC method proposed here is essentially identical to that previously described 
in Pauw et al. (2013). Differences can be summarised as: 1) the implementation of 
a different elementary scatterer: an oriented ellipsoid of revolution with semi-axes a, 
b = a and c whose axis radii will be denoted by R a and R c , respectively. 2) The use 
of the detector angle tp (c.f. in (Pauw et al, 2010)), and 3) a single MC iteration uses 
three randomly generated numbers that define R a and R c , and the in-plane rotation 
ip e \\. For computational consistency (but without a change in function), the previously 
utilised "compensation parameter p c " has been replaced by a "radius-power factor" 
w, where w = The changes in the affected MC procedural steps due to all 
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aforementioned modifications are given below. 

2.1. Scattering of an ellipsoid 

The 2D scattering of an ellipsoid of revolution onto a detector pixel characterised 
by the scattering vector length q and azimuthal angle i/> can be calculated through 
modification of the Rayleigh function for a sphere (as used in previous work (Pauw 
et al, 2010)). Here, q is defined as q = 4ir/\sm(6), where the radiation wavelength is 
denoted as A, and 26 is the scattering angle, and the ellipsoids are considered oriented 
with their axes aligned perfectly parallel to the detector plane, and its main axis c 
rotated from the detector y-axis by ip e n : 

P / p o n sm(qR c ii) - gggllCOs(ggell) m 

i'elli.q, 1, Ra, Rc) = <$ 1 „ 3 , (1) 

W-n-ell/ 

where: 

i? ell ( 7 , R a , R c ) = v /i?2 sin 2( 7 ) + ^2 cos2(7) (2 ) 

and: 

7 = tp - 1p e ii (3) 
The volume for an ellipsoid is given by: 

y cll = ^krIr c (4) 

Step 1: Preparation of the procedure 

The initial guess of the total scattering cross-section is calculated for a fixed number 
of ellipsoid contributions n s (commonly set to 500) whose radii R a and R c as well as the 
in-plane misalignment V'eii are randomly sampled from bounded uniform distributions. 
The radii sampling method is bounded by the radius range dictated by the measured 
q-range as detailed in the previous work (Pauw et al, 2013), and the orientation 
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distribution of ellipsoid axis c is bounded by 45 degrees on either side of the fibre axis. 
It is important to note that the random sampling method is not quantised, so that 
the values can therefore assume any floating-point value within the size bounds. This 
initial guess is calculated using the general equation: 

Imc(q) =-B + i E I FellM<l^k,Re,k,Rm,k) ? (§*) ' R^k^k ( 5 ) 
fc=l v ' 

where the sum runs over all ellipsoids in the finite set. B is a constant background 
term, and A is a scaling factor, which is related to the volume fraction cf> of the 
scatterers through: 



lis 



A = (pAp 2 E 4 ^pct (6) 

where Ap is the scattering contrast. The volume fraction of the scatterers <fi is defined 
as: 

, ^scatt 

<t> = — (7) 
where V nr is the irradiated sample volume, and Iscatt the total scatterer volume in 

V\rr- 

A and B are tuned to the measured scattering through optimisation using a least- 
squares residual minimisation procedure, minimising the reduced chi-squared Xr (Pedersen, 
1997), modified for 2D patterns: 



1 " 



Xr N - M j=i i- °7^Ph J w 

where N denotes the number of data points, and I meas (q,ip) and iMcil-,^) the mea- 
sured and calculated model scattering cross-section at the position defined by (q,ip), 
respectively. a(q, is the the estimated error on each measured data point, whose 
estimation method is detailed in paragraph 3.4. M is the number of degrees of free- 
dom in the fitting model, but is unfortunately ill-defined in an MC model and thus 
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set equal to two for the scaling parameter A and background contribution parameter 
B. 

Step 2: Optimisation cycle 

Step two remains the same as for the isotropic MC variant, with the exception that 
every suggested replacement ellipsoid will be defined by three new sampled random 
variables. While it is possible to suggest changing only a single parameter of the 
original ellipsoid, this may have the potential to lead to biased optimisation and is 
therefore not done. 

Step 3: Convergence and post-optimisation procedures 

After convergence has been reached, the partial volume fraction for each ellipsoid 
contribution k, <f>k is then calculated through reformulation of equation 6: 

^ = A? ( 9 ) 

where A is known through least-squares fitting from step 2. If this equation is calcu- 
lated for all n s contributions, the volume fraction of the scatterers <fi can be calculated 

as <P = E 4>k- 
k=l 

3. Experimental 

3.1. Simulation details 

A scattering pattern was simulated using lxlO 5 ellipsoids over a grid of points 
whose q and tp vectors are identical to that of the real measurement described below. 
The ellipsoids were assigned R a and R c values using a random sampling procedure 
over triangular volume- weighted distributions. The lower limit, mode and upper limit 
of the distribution for R a were 1, 1, and 20 nm, respectively. For R c , these values were 
1, 50 and 50, respectively. The orientation value for the ellipsoid axis c was randomly 
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sampled from a uniform volume-weighted distribution with lower and upper bounds 
of 80 and 100 degrees, respectively (in all scattering patterns, tp runs clockwise from 
zero at the 12 o'clock position). 

The calculated intensity was subsequently scaled to a maximum intensity of 1 x 10 6 
counts, quantised using a rounding procedure, and baseline noise of or 1 count was 
added to all calculated pixel intensity values. The data was subsequently subjected to 
a 2x2 pixel binning procedure in order to get the data fitting times to within practical 
values (typically on the order of hours for 50 repetitions). Data below q < 0.12 nm _1 
is masked to resemble the real collected data, and data above q > 3nm _1 is also 
masked due to the unequal coverage of the angular region above this value (i.e. beyond 
q rs 311m" 1 only the corner sections of the detector contain data, no longer a full 
circle) . 

3.2. Data collection details 

SAXS measurements were performed on a laboratory instrument based around a 
rotating anode generator with a molybdenum target. Mofc a emitted radiation was 
selected and focused using an Osmic confocal mirror. The incident beam was further 
collimated using three pinholes of 0.3, 0.2 and 0.45 mm in diameter, with the sample 
positioned immediately after the third pinhole. The distance from the optics to pinhole 
#1, #2 and ^3 is 0.115, 0.582 and 0.86 m, respectively. Assuming that the second 
pinhole is the beam-defining pinhole, the transverse coherence length (and therefore 
the maximum size of precipitates contributing to the scattering signal) is on the order 
of 200 nm (Veen & Pfeiffer, 2004). Nominal beam diameter at the sample position is 
estimated to be approximately 200 /im. The flightpath is a continuous vacuum from 
optics to detector. 

For detection of the scattered radiation, a Dectris Pilatus 100k detector was used 
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in combination with a translation stage. For every measurement, three images are col- 
lected and stitched together to form a single, almost square scattering image consisting 
of 487 by 497 pixels, each measuring 172 /jm by 172 p. The detector is placed 1.36 m 
from the sample, its distance verified by ruler and measurement of a silver behenate 
standard, allowing for an angular coverage of 0.1 to 4nm~ 1 in q. As the beamstop 
is mounted on thin polymer wires, there is no significant shadow from a beamstop 
holder, and thus no gap in angular coverage. 

The sample measurement was performed for 14 hours per frame, with three frames 
making up the scattering image. The background (no sample) was measured for 6| 
hours per frame. The transmission factor was measured by removing the beamstop 
and measuring the attenuated direct beam (Iq) and the attenuated beam with the 
sample in position (ii), the difference of which h/Io results in a transmission factor 
of 0.93(5). 

3.3. Sample details 

A commercially available sample of PPTA fibre (Twaron 1000) was obtained from 
Teijin Aramid B.V., with a gravimetric density of approximately 1450(20) kg/m 3 and 
an elemental composition of C7H5NO. This fibre is wound on a small rectangular 
frame of 13 by 18 mm used to standardise WAXS and SAXS measurements, as previ- 
ously reported (Pauw et at, 2010). Five such frames were packed together in order to 
maximise the scattering signal and average over a larger quantity of material, leading 
to a total thickness of aramid along the beam of 0.671 mm. As the fibres tend to carry 
moisture from the air, they were subjected to a >24h drying process in a vacuum 
desiccator prior to the measurement. 
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3.4- Data correction 

In order to achieve as high as possible relative data accuracy, the collected data is 
corrected for natural background radiation, transmission, sample thickness, measure- 
ment time, primary beam flux, parasitic background, polarisation, detector solid angle 
coverage (angular dilation) and sample self-absorption using in-house developed data 
reduction software written in the Python language. Deadtime correction is unneces- 
sary at the count rates encountered in this study (Kraft et al, 2009), and dead-pixel 
and flatfield corrections are performed by the detector acquisition software prior to 
ingestion by the data reduction software (Eikenberry et al, 2003). 

The collected intensity is binned in blocks of 2 by 2 pixels to improve calculation 
speed. Data below q < 0.12 nm" 1 and above q > 3nm _1 is masked (ibid.). 

The data is furthermore scaled to absolute units using a calibrated glassy carbon 
standard (Zhang et al, 2010). Statistical uncertainty on the datapoints is set to the 
maximum of either: 1) The propagated Poisson counting statistics, 2) 1% of the inten- 
sity in the histogram bin (the realistically minimum achievable uncertainty given the 
aforementioned corrections (Hura et al, 2000)) or 3) the standard deviation of the 
intensity values of the pixels in the bin. As the parasitic scattering is minimal due 
to complete beam path evacuation, improved statistics have been obtained by skew- 
ing the time division between sample and background measurement to favour sample 
measurement time (Pauw & Tardif, 2012). Furthermore, the minimal parasitic scat- 
tering improves the data accuracy by its reduced sensitivity to inaccuracies in the 
transmission factor determinations during background subtraction. 

3.5. MC method settings 

The number of contributions n s has been set to 500, and w has been set to 0.5. 
The scattering contrast is assumed to be the contrast between the aforementioned 
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aramid and vacuum, resulting in a scattering contrast 5p 2 of 1.64x 10 30 m~ 4 . The 
radius bounds for R a and R c for the random number generator are set to 0.7 < R < 
70 nm, the upper limit is set high as it is related to the actual beam width rather 
than the minimum measured q 2 . 50 independent converged solutions are calculated 
for each fit to obtain uncertainty estimates on the resulting distributions. 

4. Results and discussion 

The simulated data retrieval shows very good agreement of the calculated intensity 
with the simulated intensity as shown in Figure 1, fitting on average to within the data 
uncertainty (i.e. a Xr < !)• In addition, the distributions are reasonably well retrieved, 
though some of the sharper details of the distributions used for the simulation have 
been lost. This can most likely be ascribed to the interplay of the somewhat similar 
effects of aspect ratio and orientation distribution for ellipsoidal scatterers, though a 
similar loss of sharp details has also been observed in the isotropic retrieval studies, 
especially for broad distributions (c.f. SI of (Pauw et al, 2013)). As the intensity is 
well described, however, the result is considered to be a reasonable estimate of the 
actual structural parameters. 

2 Imagine, for example, a high aspect ratio prolate ellipsoidal scatterer. Its scattering pattern will 
resemble a streak in the detector plane, the width of which is (at minimum) the beam width. Thus, 
the maximum R is dictated by the beam spot width rather than the minimum measured q. 



IUCr macros version 2.1^1: 2007/05/15 



11 



Simulated vs. Fitted intensity 



Inyaruliniic scale) 



-3-2-10 1 2 



R radius histogram 



R radius histogram 



Orientation distrihution 




R jnmi 



R ( jnm) 



60 80 100 120 



Fig. 1. Leftmost: 2D Monte Carlo fit (quadrants with rounded corners) to simulated 
data, with the (simulated) "fibre" axis horizontal, and intensity shown on logarith- 
mic scale. Middle left: MC retrieved size histogram of the a-axis radius distribution 
originally simulated with lower limit, mode and upper limit of 1, 1 and 20 nm, 
respectively. Middle right: MC retrieved size histogram of the c-axis radius distri- 
bution, originally simulated with lower limit, mode and upper limit of 1, 50 and 50 
nm, respectively. Rightmost: MC retrieved c-axis orientation distribution retrieval 
from simulated data, originally simulated using a uniform distribution with lower 
limit of 80 and upper limit of 100 degrees (90 degrees is horizontal). Distributions 
shown as volume- weighted distributions. 



The Monte Carlo procedure as applied to analyse the scattering pattern effected by 
the pore structure of a PPTA fibre is shown in Figure 2. The MC method can fit the 
data on average to within the data uncertainties (i.e. with a Xr < !)• Fhe retrieved 
size distribution shows few features in the c-axis size distribution, which may indicate 
that the maximum of that distribution lie well beyond the limits of the instrument. 
The a-axis distribution shows two clear features, one distribution at around 2 nm, and 
a second one at 20 nm. This result would indicate that there is a bimodal porosity 
in the Twaron 1000 fibre. Note that the real fibre-symmetric three-dimensional void 
orientation function is approximated here using a two-dimensional c-axis orientation 
distribution, similar to the treatment described in Pauw et al. (2010). 
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Fig. 2. Leftmost: Measured data versus MC result (quadrants with rounded corners) of 
Twaron 1000, with the fibre axis horizontal and the intensity on a logarithmic scale 
in absolute units of reciprocal meters. Middle left: a-axis radius distribution. Middle 
right: c-axis radius distribution. Rightmost: c-axis alignment (fibre axis horizontally 
at 90 degrees). Error bars indicate ± 1 standard deviation from the mean calculated 
from the result of 50 MC optimisations. Distributions shown as volume-weighted 
distributions. 



In the prior evaluation of the isotropic data retrieval using the MC method, it was 
evident that shape and size distribution cannot both be uniquely identified in a scat- 
tering pattern from polydisperse systems. A similar limitation exists for anisotropic 
scattering patterns, despite the availability of additional dimensionality in the data. 
To illustrate this, two solutions were calculated for the Twaron 1000 fibre, one con- 
sisting of oblate ellipsoids, the second consisting of prolate ellipsoids 3 . Both solutions 
converge on average to within the uncertainty of the data, and are shown in Figure 3. 



3 The deciding factor for which shape was preferentially obtained is the bounds of the c-axis alignment. 
If the c-axes are allowed to vary within ±45° from the direction perpendicular to the fibre axis (and 
parallel to the detector plane), oblate ellipsoids are preferentially obtained. If the c-axis axes are 
allowed to vary within ±45° parallel to the fibre axis, however, prolate ellipsoids are obtained. 



IUCr macros version 2.1,81: 2007/05/15 



13 



Twaron fitted to oblate ellipsoids Twaron fitted to prolate ellipsoids 
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Fig. 3. Twaron 1000 measurements compared to the MC results (quadrants with 
rounded corners) for oblate ellipsoids (Left, with the c-axis within 45 degrees of the 
vertical axis) and prolate ellipsoids (Right, with the c-axis within 45 degrees of the 
horizontal axis). Intensity shown on logarithmic scale. 

These results stress the need for supplementary information when finding a solution 
to a small-angle scattering problem. In this case knowledge on the general scatterer 
shape (oblate or prolate, ellipsoids or cylinders) is needed from either microscopic 
investigations or considerations on the production process. Given the production pro- 
cess of the fibres, the existence of prolate ellipsoidal voids in the fibres is the obvious 
candidate. 

The thus retrieved data particularly lend itself to further scrutiny, as it contains 
the interrelations of both dimensions (and thus the aspect ratio) alongside the ori- 
entational information. Depending on the requirements of the researcher, therefore, 
the ellipsoid parameters resulting from the MC method can be plotted against one 
another in a variety of ways. Amongst others, this allows for the study of the corre- 
lation between the aspect ratio and particle volume, or radius distribution and ori- 
entation angle. More simple parameters can also be obtained, such as the radius of 
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gyration, the modes of the orientation distribution, and the total volume fraction of 
the phases. As such, we believe this is a very powerful and promising method of 2D 
data analysis. 

5. Upcoming improvements 

In contrast to the previously published isotropic MC method, this 2D variant can 
shamelessly be called "very slow". A single MC optimisation can take up to several 
hours on a single core of a 1.8 GHz intel i7 (as found in the 2011 Macbook Air), 
a similar timescale as the laboratory measurement itself. This appears mostly due 
to the size of the dataset (on the order of a quarter million datapoints, reduced to 
about 60000 upon 2x2 binning) and the need to repeatedly evaluate the ellipsoid 
scattering function as well as the chi-squared function with these points. One upcoming 
improvement therefore is the implementation of faster calculation functions for the 
ellipsoid scattering pattern, similar to those recently suggested by Littrell (2012). 
Furthermore, multicore processing will be implemented to leverage more computing 
power. 

Secondly, the issues arising from the interplay between the similar effects of aspect 
ratio and orientation distribution may be further suppressed when cylindrical shapes 
are used instead. Cylindrical shapes have scattering patterns which are more direction- 
dependent than their ellipsoidal counterparts, and contributions arising from their 
aspect ratio should therefore be better separable from the orientation distribution 
effects. Alternatively, a Bayesian approach could be considered in which prior infor- 
mation from other techniques is used to arrive at a inter-technically consistent solution. 
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6. Conclusions 

This report shows the applicability of MC methods to anisotropic scattering patterns, 
which can be used to extract all available information from a scattering pattern. A 
wide range of morphological questions can be answered through careful interpretation 
of the MC result. While the initial results are encouraging, improvements abound, in 
particular on improving the speed and accuracy of the method. 

The MC code as well as the measured datasets are available for inspection, improve- 
ments and application under a Creative-Commons Attribution Share-alike license and 
the latest copy will be freely supplied by the author upon request. 
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